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Abstract 

Analytical and numerical modelling of the nonlinear interaction of shear wave with a 
frictional interface is presented. The system studied is composed of two homogeneous and 
isotropic elastic solids, brought into frictional contact by remote normal compression. A 
shear wave, either time harmonic or a narrow band pulse, is incident normal to the interface 
and propagates through the contact. Two friction laws are considered and their influence 
on interface behavior is investigated : Coulomb's law with a constant friction coefficient and 
a slip-weakening friction law which involves static and dynamic friction coefficients. The 
relationship between the nonlinear harmonics and the dissipated energy and their depen- 
dence on the contact dynamics (friction law, sliding and tangential stress) and on the normal 
contact stress are examined in detail. The analytical and numerical results indicate universal 
type laws for the amplitude of the higher harmonics and for the dissipated energy properly 
non-dimensionalized in terms of the pre-stress, the friction coefficient and the incident am- 
plitude. The results suggest that measurements of higher harmonics can be used to quantify 
friction and dissipation effects of a sliding interface. 

1 Introduction 

The detection of contacting interfaces which can slide or open under applied load is a difficult 
task but one that is relevant in many technologies. Examples include closed cracks or defects, 
glued bondings, and partially delaminated material interfaces. Ultrasonic methods based on 
linear wave scattering are efficient for detecting defects and characterizing material elasticity, 
but are less sensitive to closed cracks and contacting interfaces. Using the nonlinear behavior of 
these defects, nonlinear ultrasonic techniques such as nonlinear resonance [1], sub- and higher- 
harmonic generation [2] , and frequency-modulation [3] have been shown to be sensitive to micro- 
cracks or closed cracks. For an overview of nonlinear acoustics applications see [UEIE]. When an 
ultrasonic wave with large enough amplitude is incident on a contact with a frictional interface 
(e.g. closed cracks), higher harmonics appear in the frequency spectrum of transmitted and 



reflected waves. This effect, called Contact Acoustic Nonlinearity (CAN) [2J, is of increasing 
interest for characterization of closed cracks or imperfectly bonded interfaces [HE]. 

Analysis of model systems has helped to understand the nonlinear interaction of waves and 
frictional cracks. In particular, [8J analysed the partition of energy resulting from a time har- 
monic obliquely incident plane SH wave reflected and refracted by a frictional contact interface 
that slips. Both longitudinal and shear wave propagation through a rough surface were inves- 
tigated in [9J using an interface contact model based on Hertz theory, again for time harmonic 
incidence. The partial contact model was subsequently applied to model scattering from surface 
breaking cracks [10J, and numerical simulations indicated efficient production of second har- 
monics. Time domain studies have concentrated on numerical implementations, such as BEM 
formulation of SH slip motion on an arbitrary interface [llj. Using a generalization of this 
method to include in-plane motion, it was shown in jT2j [13] that the amplitudes of the higher 
harmonics of the scattered far-fields can be useful in determining both the pre-stress and the 
frictional coefficient. Measurements of second harmonic generation for normal incidence of longi- 
tudinal waves on a contacting interface between aluminum blocks have been reported in [HI [15]. 
These experiments indicate that the amplitude of the second harmonic decreases rapidly with 
applied normal contact pressure initially, and then falls off in magnitude at a lesser rate. These 
findings are in agreement with experimental measurements on contacting adhesive bonds |16| . 
A simple but instructive model was proposed in [17J for the nonlinear interaction of an SH wave 
normally incident on a frictional contact interface. In contrast to the other studies mentioned, 
this friction model predicts that only odd harmonics are generated for time harmonic incident 
wave motion. 

The objective here is to investigate the effect of friction on shear wave scattering from a 
interface. Physically the friction leads to nonlinearity due to switching between the sticking 
and sliding states of the interfacial contact. Analytical and numerical methods are combined to 
understand the role of energy loss and the generation of higher harmonics and their dependence 
on the system parameters, such as applied normal load, incident amplitude, and the friction 
coefficient. We consider a shear wave incident normally on a flat frictional interface. The friction 
behavior of the interface is described by two friction laws: Coulomb's law with a constant friction 
coefficient and a slip- weakening friction law including static/dynamic friction coefficients. The 
outline of the paper is as follows. The analysis in Section [2] develops the basic model of [17] to 
consider both time harmonic and transient pulses. Particular attention is given to estimating the 
energy lost in frictional sliding, and how it effects the scattered waves. A closed form expression is 
obtained for the energy dissipation under time harmonic incidence. The time domain numerical 
model is presented in Section |3j We note that the numerical simulation contains more of the 
actual physical effects encountered in practice, such as finite contact interface and cylindrical 
wave spreading. Simulation results are discussed in Section [4] in light of the simpler analytical 
model, and several conclusions are drawn concerning the relation between the total energy 
dissipated and the amplitudes of the harmonics in terms of non-dimensional system parameters. 

2 Nonlinear analysis of shear wave interaction 
2.1 Setup and energy dissipation 

We consider uni-dimensional SH motion in a solid with shear modulus G, density p, impedance 
Z = (Gp) 1 / 2 and wave speed c = {G/p) 1 ' 2 . A plane SH wave is normally incident from y > 



2 



on the interface y = 0, such that the total shear stress a xy = a T comprises the incident and 
reflected waves above the interface, and a transmitted wave below: 

{a inc + a re f, y > 0, cr inc = f(t + |), 

' where a ref = g(t - ^), (1) 
Gt ^ V< ^ ci tra = h{t+l). 

Instantaneous particle velocity in the x direction is v(t,y) and equal to Z~ 1 (ai nc — a re f) for 
y > and Z~ x o tra for y < 0. Shear stress is continuous across the interface, regardless of the 
state of contact, with stress <j t (£, 0+) = cr T (t, 0—) = cr T0 (t), where the latter follows from as 

a T0 (t) = f(t)+g(t) = h(t). (2) 
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Figure 1: (a) The filled part depicts the transmitted wave for an incident time harmonic wave of 
period T and amplitude twice the critical value a TC = 0.5. The remaining part of the wave form 
between h and / is the negative of the reflected wave g(t), by Eq. (b) A narrow band pulse 
for the incident stress function f(t). The central frequency is uo — 2tt/T with T = §, a TC = 0.3. 

When the interface is in contact the wave is completely transmitted with zero reflection: 
g(t) = 0, h(t) = /(£), and there is no energy dissipation. When slipping occurs, the velocity 
jump across the interface follows from ([!])- ^ as 

[v(t)\ = v(t, 0+) - v(t, 0-) = -2Z- X g{t). (3) 

Slipping incurs irreversible loss of mechanical energy into heat, which can be quantified in two 
equivalent ways. First, the work done per unit interface length by the slippage over a time 
increment At is (force x distance =) a T0 [v]At. The total amount of energy dissipated (per unit 
interface length) then follows from eqs. ^ and ^ as 

D = [ <T T0 (t)[v(t)]dt = -2Z' 1 [ g(t)h(t)dt. (4) 
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(a) (b) 

Figure 2: (a) Magnitudes of the third and fifth harmonic amplitudes of the reflected or transmit- 
ted wave. The third curve is D/1Q where D is the nondimensional energy dissipation function 



of (10). The abscissa is £ = ii\a N \/a TM = a TC /a TM , where a TM is the maximum stress of the 
incident shear wave. The maximum value of As is at £ = 0.5, and maximum of D occurs at 
£ = 0.3942. (b) Reflected and transmitted first harmonic amplitudes. 

Alternatively, the rate of flux of incident energy per unit length (power flow) is (stress x velocity 
=) Z~ x o 2 nc . The fluxes of the reflected and transmitted waves are Z~ x a 2 e ^ and Z~ 1 a 2 ra , re- 
spectively. The loss in energy (per unit length and time) is therefore Z~ 1 (a 2 nc — cr 2 e j — of m ). 
The three waves are assumed to be plane waves, and hence the fluxes are independent of y, with 
y > for the incident and reflected waves, and y < for the transmitted wave. Taking y = ±0 
and using to express the stresses, the total energy loss per unit length can be written 

wave energy lost = Z' 1 J (f 2 (t) - g 2 (i) - h 2 (t)) dt. (5) 

Substituting from eq. ^ shows that the lost wave energy is the same as the dissipated interfacial 
energy of eq. Q. 

2.2 Coulomb friction 

Before the incident wave strikes the interface it is assumed to be in a state of constant normal 
stress a N < (compression). Slipping commences at the first instant that 

|<j T0 (t)| = a TC , (6) 
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where the critical value a TC depends on the friction model. For simple Coulomb friction: a TC — 
n\cr N \, where fi > is the frictional coefficient. Slipping continues for as long as \f(t)\ > a TC 
and \f(t) + g(t)\ = a TC while it occurs. By considering the two possibilities: f(t) > a TC , 
f(t) < —cr TC , it follows that the reflected wave during slipping is 

9(t) = -(1/(01 - cr TC )sgn(/(t)) H(\f(t)\ - a TC ), (7) 

where H is the Heaviside step function. The transmitted wave amplitude h is given by eq. ^ 
and the dissipated energy is therefore 



D = 2Z- 1 a rc J (\f(t)\-a TC )dt. (8) 

t:\f(t)\>a TC 

The slip condition Q is equivalent to the model proposed in [T7]. It predicts that the reflected 
and transmitted waves are "clipped" versions of the incident pulse at the positive and negative 
values equal in magnitude to the critical yield stress, as illustrated in Figure [T] The simple 
expression Q for the dissipated energy is new and will be central to subsequent developments. 

As an example, consider the time harmonic incident wave f(t) = a TM smujt, a TM > a TC . 
During a single period, < t < T = 2tt/uj the slip condition \f(t)\ > a TC prevails for t E 
{tc, \ — tc) and t G (|+ tc-, T — tc) where uotc = 0c = sin _1 (o rc /a rM ) E (0,7r/2), with the 
same amount of energy lost in each sub-interval. It follows that the reflected stress function is 

-g^Q) =a TM (sme-sm6c)H(e-e c ) for < 9 < it/2. (9) 

The energy dissipated in one cycle then follows from Eq. Q as 

D = ^fD(Z), £=— , ^(O = 8C[(l-C 2 ) 1/2 -^(?- sin_1 0] 5 (10) 
cjZ a TM ' 2 

where the nondimensional function D is shown in Figure [2J Nonlinear effects occur only for values 
of the non-dimensional stress parameter £ less than unity. If £ > 1 no energy is dissipated as the 
wave transmits unaltered. We will also find the scaling factor ujZa~^ critical for interpreting 



the numerical simulations (see £3.3). 



2.3 Generation of higher harmonics 

The reflected and transmitted stress waves together sum to equal the wave shape of the incident 
wave. This means that the the higher harmonic amplitudes of the reflected and transmitted 
waves equal the content of the incident wave. If the latter is a single frequency time harmonic 
wave, then higher harmonic amplitudes of the reflected and transmitted waves are equal in 
magnitude and of opposite sign. We consider this case first. 



2.3.1 Time harmonic wave motion 

The incident wave is f(t) = a TM sincjt, a TM > a TC . The clipped nature of the reflected wave 
implies that its Fourier series contains only odd harmonics [17] . i.e. 

00 2 2/ 

9(t) = <y TM V A 2n _ l sin((2n - A 1 &) = 1-- sin" 1 £ - - £ 2 ) 1/2 , 

Z / TV 

71=1 

As(0 = -^(l-ef\ A(0 = -7^(3- 8^(1 -£ 2 ) 3/2 > etc., (11) 



5 



where £ is defined in (10) . The first harmonic amplitudes and the magnitudes of the third and 



fifth harmonics are plotted in Figure [2j Note that by definition the non-dimensional amplitudes 
are normalized such that A\ + B\ = 1 where B\ is the amplitude of the first harmonic of the 
transmitted wave. 
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Figure 3: The dissipation and the second, third and fifth harmonic amplitudes as a function of 
£ = Vtc/^tm wnere &TM i s ^ ne amplitude of the incident wave in Figure [TJ). 



2.3.2 A narrow band pulse 

A narrow band pulse is shown in Figure [TJ). Note that the pulse is antisymmetric about the 
center t = because it represents the incident shear stress, which is proportional to the velocity. 
The incident displacement amplitude is symmetric about t — 0. Explicit expressions for the 
dissipation and for the harmonic content are not available, but are easily found by numerical 
quadrature. Thus, the dissipation is computed from Q, and the harmonic amplitudes of the 
reflected wave are given by analogy with (11)2 as 



A r . 



UJ 



7T 



a TM / g(t)sm(nujt)dt. 



(12) 



Note that all of these quantities (D, A n ) are proportional to Fourier components of the reflected 
wave, either at zero frequency for D, or finite values for the harmonics. Figure [3] shows the dissi- 
pation (normalized by the scaling factor uZa~^ of Eq. (jToj)) and the magnitudes of the several 
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harmonics. Comparing Figures [2^i and [3] we note that the harmonics A3 and A$ for the pulse 
are larger in magnitude as compared with the time harmonic case, whereas the nondimensional 
dissipation D is smaller. The overall magnitudes are still comparable with those in the time 
harmonic case because, even though the integral is over all t rather than a single period, the 
majority of the wave clipping and dissipation occurs within one cycle. Note that A2 7^ but it 
is small in comparison with A3. 

We next consider a more realistic model using numerical simulation. 

3 Numerical study for non-plane wave incidence 
3.1 System studied 



Imposed displacement 
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Figure 4: (a) The system considered in the numerical model, (b) Snapshot of the velocity 
magnitude illustrating the non-plane wave incidence generated from the upper surface of pW. 

The system modeled is shown in Figure [4j It is composed of two identical isotropic elastic 
media and P^ 2 ) separated by an interface governed by unilateral contact with a Coulomb 
friction law. The material used in the simulations is aluminium with the Young modulus, 
Poisson's ratio and density equal to 69GPa, 0.33 and 2700kg/m 3 respectively. The problem is 
considered in terms of a material (Lagrangian) description. 

In the first step of the simulation a normal contact stress between the two solids is applied 
by numerically imposing appropriate boundary conditions on the external surfaces. The lower 
surface of the body P^ 2 ) is fixed (u^ = 0) and the upper surface of the body P^ is shifted by 
the appropriate amount Vq in the y direction (u^ — Vq). The normal contact stress profile is 
shown in Figure [5^. The normal contact stress a N is quasi-constant on a part of interface (in 
Figure^, a°=-0.3MPa). 

In the second step, a non-planar shear wave is generated by imposing a dynamic displacement 
along the x axis on the upper surface of P^ 1 ). During this stage of the simulation the lower surface 
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Figure 5: (a) Normal contact stress along contact interface after the first step of simulations 
for \a®\ = 0.3MPa. (b) The imposed displacement on the upper surface = W(x)U(t). 

of body P^ 2 ) is fixed (u^ = 0) and the displacement imposed on the upper surface of P^ 1 ) is 
expressed as a product of a weight function W(x) and a function of time U(t) shown in Figure 
[5J3. The interface is supposed to remain in contact during the transmission/reflection of the 
shear wave. This assumption has been verified for all simulations. 

3.2 Numerical model 

PLAST2, an explicit dynamic finite element code in 2D, is used to simulate the behaviour of the 
system (pW and P^ 2 )) during frictional contact. PLAST2 is designed for large deformations and 
non-linear material behaviour |18j . It uses a forward Lagrange multiplier method for the contact 
between deformable bodies. For this dynamic study, the formulation is discretized spatially using 
the finite element method and discretized temporally by the 02 method. The contact algorithm 
uses slave nodes (situated on the contact surface of pW) and target surfaces (on the contact 
surface of P^ 2 )) described by a four node quadrilateral element with 2x2 Gauss quadrature rule. 
The elementary target segments are described by two nodes and approximated by bicubic splines 

The forward Lagrange multiplier method is formulated for the equations of motion at time 
t % — iAt with the displacement condition imposed on the slave node at time £ z+1 : 

MLP + CtT + KIP + G i+lT A* = F, G i+1 U i+1 < 0, (13) 

where M, C and K are symmetric and positive definite matrices of mass, Rayleigh's proportional 
damping (C = diM + cfoK) and stiffness of the system, respectively. U, U, U are respectively 
the vectors of nodal displacements, nodal velocities and nodal accelerations. F is the vector 
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of nodal external forces. The vector A = [\ni^t] T comprises normal (N) and tangential (T) 
forces at contact points C m . G T = [G^, G^] is the global matrix of the displacement conditions 
ensuring non-penetration and the contact law of the bodies in contact. 
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Figure 6: Tangential stress as a function of slip for Coulomb friction law (a) and for slip- 
weakening friction law (b) . 



The equations of motion (13) are discretized in time using an explicit Newmark scheme. The 



vectors are expressed at each time step using a time scheme of type /?2, (P2 £ [0.5, 1[): 



^(U^-lP-AilP), 



^(U^ 1 + At(l - /3 2 )U i " 1 + 2 -£(XJ i+1 



IP)). 



(14) 



The displacements *IP of the nodes situated on the contact surface (P^ and P^ 2 )) are first 
computed with X 1 equal to zero. For simplicity /?2 is fixed at 0.5 and the nodal displacements 
*IP +1 at time t l+1 are obtained according to 



*LP +1 = At 2 M _1 (F - KIT) + 2U* - U <_1 . 



(15) 



A constraint matrix G t+1 is formulated for the slave nodes if they have penetrated through a 
target segment. Calculation of the contact forces A* and of the correct nodal displacements U* +1 
at time t %+1 are then performed: 



'\ i = (Af 2 G i+1 M _1 G i+ir ) 1 G i+1 *U i+1 , 
U*+i = *ui+i _ At 2 M~ 1 G* +lT A i . 



(16) 



Equations (16) are solved using the Gauss-Seidel method. The contact conditions solved during 



each internal iteration of this method are expressed as follows : 



& N < 0, (contact if cr N < and separation if a N = 0) 

i ,^ I I a I v t=° H \cr T \ < fi a \a N \ (stick), 

W T \ <MakJ, and <^ (17) 

I (JrpVrp < II \(T T \ = Ha\V N \ (Slip), 
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Figure 7: Propagation of shear wave in (a) the body P^ 1 ) and (b) in P^ 2 ), both for x = 0. The 
waterfall plots show the x-velocity versus time for three different friction coefficients (/i= 0, 0.4 
and 1). FFT of x- velocity of (c) reflected and (d) transmitted waves for three different friction 
coefficients (fi= 0, 0.4 and 1). The FFT processing of the reflected and transmitted wave are 
respectively done at point A (x = 0,y — 20mm) and point B (x = 0,y — —5mm). 
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where v T is the relative tangential velocity for a slave node related to the surface, a N and a T 
are normal and tangential stress respectively, subscript N and T referring to the normal and 
tangential vectors defining the contact. fi a is an apparent friction coefficient. Two friction laws 
are considered in the simulations : a standard Coulomb friction law with a constant friction 
coefficient (/x a =/x, see Figure [6^1) and a slip- weakening (SW) friction law (Figure [6)3) for which 
[la is defined in [20j [21] as: 

*.M- h as) 

I /id, Uj, > a, 

where jig and \±b are static and dynamic friction coefficient, is the slip and d the critical slip 
distance. Once the dynamic regime is reached, the friction corresponds to a classical Coulomb 
friction law with a dynamic friction coefficient fio (< fJ>s). This physically motivated model of 
slip weakening friction (i.e. the decrease of the friction force with slip), which was first introduced 
in the geophysical context [22], can be associated to a phenomenon of adhesion caused by long 
time stationary contact [23]. For more details on rate and state- variable friction law (including 
the slip- weakening one) see for example |24| . 



3.3 Results 

3.3.1 Case of Coulomb friction law 

Figure [7] shows the velocity on the x-axis versus time and y-position (depth) for three different 
friction coefficients. At time t = 0, the shear wave is generated on the upper face of P^^see 
Fig. [4]). The incident wave then propagates through pW and the contact interface gives rise to 
reflected and transmitted waves. 

It can be seen from Figure [7] that the magnitude of the transmitted wave is increasing with 
the friction coefficient. Consequently, the magnitude of the reflected wave is decreasing. Points A 
(x=0,y =20mm) and B (x=0,y=-5mm) are chosen for Fast Fourier Transform (FFT) processing 
of the reflected and transmitted waves, respectively. 

Figure^ and[7]l shows the FFT of velocity along the x-axis for the reflected and transmitted 
wave for /x=0, /x=0.4 and fi=l. Nonlinear harmonics should appear in the spectrum as soon 
as slipping zones occur at the contact interface. For friction coefficient equal to or 1, no 
harmonics are observed in the reflected and transmitted wave spectrum. For the case fi— 0, 
the interface is in a sliding state during transmission of shear waves but the system is perfectly 
linear under the previously made assumption of no separation. When the friction coefficient 
is equal to 1, the contact law is nonlinear, but the friction coefficient is too large to permit 
enough sliding such that would produce a noticeable effect on the transmitted/reflected wave 
spectrum. In other words, during transmission/reflection of the shear wave the contact nodes 
remain almost always in the sticking state of the Coulomb friction law. For friction coefficient 
equal to 0.4, the third harmonic (1.5MHz) is observed in the reflected and transmitted wave 
spectra due to the nonlinearity of the friction law and the occurrence of sliding zones during the 
transmission/reflection stage. 

As expected, the magnitude of the wave decreases along the propagation direction because of 
diffraction effects. In order to compare the numerical results with the analytical ones presented 
m 

g the non-dimensional parameter £ = fi\a N \/a TM used in the analytical model (where a TM 
denotes the maximal stress magnitude for the plane wave) needs to be modified. Let £ = 1 
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correspond to the case in which the wave magnitude is not large enough to initiate slip at 
the interface. From this definition a critical friction coefficient for slip, \i C riti is determined 
numerically, and the modified £ num is defined as follows: 



Cnum 



O num. 



with G, 



num — l^crit]®. 



I 



(19) 



The harmonic magnitudes are extracted from the FFT of the computed velocity (Figure [7^ 
and d). Energy dissipated during transmission of the shear wave D(x) is numerically calculated 
along the interface from its definition Qi. Since the incident wave is non-plane, D(x) is not 
constant along the interface, and accordingly a mean value D on the zone of constant normal 
stress (Figure [5^i) is chosen for the results. Non-dimensional measures of the harmonic magni- 
tudes and energy dissipated per unit cycle during transmission of the shear wave are defined 
respectively as 

- * - » Z ~ (20) 



Ai 



J-r? 



D 



hoi 



-D, 



where A¥ is the magnitude of the first harmonic of the transmitted wave velocity for £ = 1 
at the probe point B and the factor 5 is the number of cycles in the incident signal (see Figure 
Kb for U{t)). 
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Figure 8: (a) Third and fifth harmonics A3, A§ and dissipated energy D for \a®\ = 0.3MPa 

( ) and for = 0.6MPa (x). (b) Tangential stress and relative displacement on interface 

(x=0) versus time for three values of ^ num . The results are presented for the case \a® |=0.3MPa. 

Figure [8^1 displays the evolution of the non-dimensional third and fifth harmonic magnitudes 
and non-dimensional dissipated energy versus ^ nU m fc> r two different values of initial normal 
stress: 0.3MPa and 0.6MPa. It can be seen that the initial normal contact stress has no influence 
on the non-dimensional third, fifth harmonics and dissipated energy. Comparing Figures [8^1 and 
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[3| shows a good qualitative agreement between the numerical and analytical results for D and A{. 
The fact that the values obtained numerically are about half of the analytical ones is because 
the sliding is not uniform for a cylindrical shape of wave front used in the numerical model, 
while the analytical model assumes a plane wave. 

The appearance of odd harmonics in the spectrum is due to distortion of the incident wave 
during the transmission/reflection as the interface partially slides under friction. The evolution 
of the odd harmonics, directly related to this distortion, depends on two competing effects. 
First, for a fixed value of the sliding distance, more energy is dissipated and the wave is more 
highly distorted as the friction coefficient increases. On the other hand, since the incident wave 
is the same for the different simulations, the cumulative sliding distance is actually diminished 
as the friction coefficient increases, as shown in Figure ^jp. These two antagonist aspects of 
nonlinearity of a frictional interface explain the existence of a distinguished value of the normal 
stress that maximizes energy dissipation: ^ num ~ 0.4 in Figure^. 

The numerical simulations allow us to explore dependence on parameter values describing 
the contact interface (displacement, normal and tangential stresses, sliding distance,...). Figure 
^jp shows the tangential stress and displacement of the contact node situated on the interface 
at x — mm versus time. The displacement corresponds to the tangential displacement of the 
slave node (of pW) in a local frame of P^ 2 ) whose origin is situated at the initial contact position 
of this node (of pW) on the P^ 2 ) interface before the arrival of the incident shear wave. As the 
friction coefficient increases the relative displacement decreases, because the resistance to sliding 
is increasing. For £ num =0.1, the sliding of the slave node is important but the sinusoidal form 
is weakly affected, whereas for ^ nU m =0.8 the sliding of the slave node is weak, but the evolution 
of the displacement of the contact node is strongly distorted from the initial waveform. 

3.3.2 Case of slip-weakening (SW) friction law 

In this simulation the numerical model takes into account a slip-weakening friction law that 
includes [201 EI] the combined effects of static and dynamic friction, with coefficients /is and 
Hd respectively. A third parameter d describes the slope of the decreasing friction coefficient, 
see Figure [6)3. For all simulations, static friction coefficient and incident shear wave magnitude 
remain the same ( /is=0.3 and a TM =0.33MPa). The simulations were performed for different 
values of the initial normal stress \cr® | (from 0.15 MPa to 1.2 MPa). Note that in these conditions 
no slip is allowed when is greater than l.lMPa. 

Figure [9] shows the results for two dynamic friction coefficients in comparison with those 
for the constant friction law. Figure [9^i confirms that for this model the third harmonic again 
follows the behavior of the dissipated energy. Maximum values of the third harmonic and 
dissipated energy per unit surface D are decreasing when the dynamic friction coefficient is 
decreasing. The maximum of both is obtained for a larger value of |<t^|. 

Figure [9|3 shows the tangential stress a T and the tangential slip of the contact node 
situated at x = on the interface versus time. The interface is sticking because |<r T | is less 
than /is I o^J I . When |<j t | is equal to /is|a^|, the system begin to slide. |<r T | remains the same for 
constant Coulomb friction (/i = /is) while, for the SW law, \a T \ = lidW®\ after the slip distance 
becomes greater than d. Hereafter the interface behavior corresponds to a classical Coulomb 
friction law with /i = /ip. As expected, since the dynamic friction coefficient is less than the 
static one, the relative displacement at the interface is greater than for the constant friction 
coefficient case (/i = /is). Also, the relative displacement increases when the dynamic friction 
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Figure 9: (a) Evolution of the third harmonic A\ and of the dissipated energy per unit surface 
D as a function of \a®\ for Coulomb and SW laws for different friction coefficients ((is — 0.3, 
d = 10 _2 /xm). (b) Tangential stress and displacement on the interface (x — 0) versus time for 
0.3MPa. 



the case |cr° 



coefficient decreases. The resistance to motion due to friction is less important and so the 
dissipated energy is reduced. This explains the offset of the position of the maximum dissipated 
energy D (Figure^). The form of the curve A3 ( | cr^ | ) is directly related to the interface friction 
law in effect. 



Figure 10 highlights an important effect of the parameter d on the tangential stress and slip 
at the interface, and consequently on dissipated energy and third harmonic magnitude. For a low 
value of d (= 10 _3 /xm), the system behavior is similar to the constant friction law with (i — \ij^. 
The sliding is delayed due to the static friction coefficient. For low value of |o^|(< O.QMPa) 
this delay has little influence on the interface behavior. On the contrary, if > 0.6MPa, 
the delay is growing and the interface slip decreases reaching zero for |cr°| = l.lMPa. This limit 
value of I cr^l depends on the parameter d namely, it decreases as d increases. As a result, for 
d = 10 - Vm, the limit value of is very low and hence the dissipated energy and the third 
harmonic curves tends to those for constant friction coefficient /i = [is- 



4 Conclusions 



This study presents a combined numerical and analytical analysis of the propagation of a shear 
wave through a frictional interface. A simple analytical model is developed that describes the 
case of a classical Coulomb friction law and a plane wave. The model pinpoints the pivotal 
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Figure 10: (a) Evolution of the third harmonic A\ and of the dissipated energy per unit surface 
D as a function of \cr®\ for Coulomb (n — 0.3 and fi = 0.2) and SW laws for different values 
of parameter d ({is = 0.3, {id — 0.2). (b) Tangential stress and displacement on the interface 
(x=0) versus time for the case \a®\ = 0.3MPa. 



non-dimensional parameters of the problem and yields an explicit dependence of the measurable 
quantities on these parameters. Good agreement between the analytical and numerical results 
for this case is demonstrated. The numerical analysis is also applied to a more complex case 
which involves a non-planar wave and a slip- weakening (SW) friction law with static and dynamic 
coefficients {is / Hd- For the both cases, it is shown that the amplitude of the third harmonic is 
directly related to the friction-induced dissipated energy. The evolution of the third harmonic 
as well as of the dissipated energy depends on two aspects: on the cumulative sliding distance 
during the transmission of the shear wave and on the friction coefficient. When a SW friction law 
is taken into account, the maximum values of third harmonic magnitude and of the dissipated 
energy are achieved for a larger value of the applied normal stress |<t^|. This is explained by the 
fact that the static effect in the SW law delays the interface sliding in comparison with the case 
H = hd- Less sliding and {id < l^s both imply lower values of dissipated energy and of third 
harmonic. 

The main conclusions of this study are as follows : 

(i) There exists a specific value of the normal stress that maximizes energy dissipation. This 
value results from the competition of two antagonist aspects of a frictional interface: sliding 
(cumulative sliding distance) and resistance to motion (tangential stress). 

(ii) The SW friction law affects the form of the curves for the dissipated energy, D(\a®\), and 
third harmonic, A^dcr^l). The maximal value is obtained for a larger value of |<t^|. Each of 
the three parameters of the SW friction model has a different impact on interface sliding and 
dissipated energy: the dynamic friction coefficient {ijy acts on the steady state sliding, the static 
friction coefficient {is delays the sliding in comparison with the constant friction law /i = 
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and d acts directly on this delay. 

(iii) For the form of shear incident wave used in this study, the third harmonic evolution follows 
the variation of the dissipated energy at the interface. Thus the measurement of transmitted 
or reflected wave signal can provide information on dissipated energy at the interface, and the 
form of the A3 ( | cr^ | ) curve can reveal the type of friction law of the interface. 
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